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ABSTRACT 

By comparing near-infrared spectra with atmospheric models, we infer the effective temperature, 
surface gravity, projected rotational velocity, and radial velocity for 21 very-low-mass stars and brown 
dwarfs. The unique sample consists of two sequences in spectral type from M6-M9, one of 5-10 Myr 
objects and one of >1 Gyr field objects. A third sequence is comprised of only ~M6 objects with 
ages ranging from <1 Myr to >1 Gyr. Spectra were obtained in the J band at medium (R~2,000) 
and high (R'^20,000) resolutions with NIRSPEC on the Keck II telescope. Synthetic spectra were 
generated from atmospheric structures calculated with the PHOENIX model atmosphere code. Using 
multi-dimensional least-squares fitting and Monte Carlo routines we determine the best-fit model 
parameters for each observed spectrum and note which spectral regions provide consistent results. We 
identify successes in the reproduction of observed features by atmospheric models, including pressure- 
broadened K I lines, and investigate deficiencies in the models, particularly missing FcH opacity, that 
will need to be addressed in order to extend our analysis to cooler objects. The precision that can be 
obtained for each parameter using medium- and high- resolution near-infrared spectra is estimated 
and the implications for future studies of very low mass stars and brown dwarfs are discussed. 
Subject headings: infrared: stars — open clusters and associations: individual (Upper Scorpius, TW 

Hydrae Association) — stars: atmospheres — stars: low-mass, brown dwarfs — 

techniques: spectroscopic 



1. INTRODUCTION 

Determining the physical properties of brown dwarfs 
and very low- mass stars is important for our understand- 
ing of a broad range of topics including star and planet 
formation, circumstellar disks, dust formation in cool at- 
mospheres, and the initial mass function. The direct 
measurement of mass and/or radius for a brown dwarf 
or very low-mass stars is possible only for certain binary 
systems, which are rare^. Determining other physical 
properties like effective temperature and surface gravity 
from bolometric luminosity estimates requires a precise 
distance measurement and assumptions about age (to 
determine effective temperature) plus mass and radius 
(to infer surface gravity from evolutionary models). Al- 
though radius is not expected to change very much for 
objects older than a few 100 Myr, reliance on evolution- 
ary models is problematic as they are poorly calibrated 
at young ages. Commonly used conversions fro m spec- 
tral t ype to effective temp erature (e.g. Golimowski et all 
120041 for field objects and Luhm an et al.i ,2003 for young 
M dwarfs) rely on model isochrones and monotonic re- 
lationship between spectral type and effective temper- 
ature, which might not accurately represent the com- 
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plicated formation and evolu tion of very-low-mass ob- 
jects fe.g. lStassun et al]|2006[) . Evolutionary models also 
require atmosphere models to provide boundary condi- 
tions. Thus, synthetic spectra from such atmosphere 
models potentially offer a more direct method of inferring 
physical properties by comparison with observed spectra. 

The NIRSPEC Brown Dwarf S pectroscopic Survey 
(BDSS, iMcLean et all 120031 l2007f) is a large sample 
of high-quality near-infrared spectra obtained with the 
Keck II 10-meter telescope. One of the goals of the BDSS 
is to provide as direct a measurement of effective tem- 
perature, surface gravity, and metallicity as is possible 
for non-eclipsing objects. In principle this can be accom- 
plished by comparing observed infrared spectra with syn- 
thetic spectra from atmosphere models, without relying 
on spectral type, distance measurements, or estimates 
of radius, mass, and age. Other authors have applied 
this promising but as yet imperfect technique to brown 
dwarfs. The primary difficulty of applying this method 
to low-mass stars and brown dwarfs is the complexity 
of the spectra to be modeled, leading to imposed limi- 
tations on resolution, wavelength range, and/or sample 
size. For example, Mohanty et al. ( 2004a, j) used narrow 
wavelength ranges at high resolution and lCushing et aTl 
(Epo s,) used a broad wavelength range at low resolution. 
Reiners et al.l ()2007D combined high-resolution spectra 
and broad wavelength coverage for a small sample of ob- 
jects. 

Many uncertainties remain in the determination of 
physical properties from the comparison of observed and 
synthetic spectra so it is important to establish the 
utility of this method using a larger sample of objects 
and a broader wavelength range. We accomplish this 
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(R~20,000) J-band spectra from the extensive BDSS 
data base^. Our goal is to show the extent to which syn- 
thetic spectra can be used to derive physical properties 
(effective temperature [T^ff], surface gravity [log(g)], 
projected rotational velocity [wsini], and radial veloc- 
ity [RV)]) rather than the measurement and analysis of 
trends in the observational data. This paper focuses on 
field (> 1 Gyr) and young (1 to <100 Myr) M dwarfs, the 
majority of which are predicted to be substellar. Future 
papers will deal with a larger sample of young M dwarfs, 
young and field L dwarfs, and finally T dwarfs. 

Section 2 describes the target selection, NIRSPEC ob- 
servations, and data reduction. Section 3 outlines the 
calculation of atmospheric structures and synthetic spec- 
tra with the PHOENIX model atmosphere code and our 
novel spectral fitting method. Section 4 presents the 
analysis of the spectral fits, including an overview of 
the results for each resolution and wavelength range and 
for the inferred physical properties. Our results are dis- 
cussed in § 5. Section 6 summarizes our conclusions and 
the implications of this work for future studies of brown 
dwarfs. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Sample Selection 

The sample of 21 objects constitutes three overlapping 
subsamples: 1) young (~5-10 Myr) objects with spec- 
tral types from M6 to M8.5, 2) field objects (> 1 Gyr) 
with spectral types from M5.5 to M9, and 3) ~M6 ob- 
jects with ages from <1 Myr to > 1 Gyr. Targets were 
selected based on pubhshed estimates of spectral type 
(M5.5 to M9), brightness (J < 14.5), and age. The 
names, spectral types, coordinates, and J-band magni- 
tudes of the sample are listed in Table [TJ Often the sam- 
ple contains more than one object of a given spectral type 
and age to look for possible differences between objects 
(effects of multiplicity, rapid rotation, etc.) and to aid in 
identifying fiducial objects with well-determined physi- 
cal properties. Ages of young objects are determined 
from membership in near-by star-forming regions and 
moving groups or confirmed companionship to known 
young stars age ~1 Myr to <100 Myr, with the excep- 
tion of 2MASS 0608-27, which has not been definitively 
associated with any other young objects. The young 
spectral type sequence contains objects with ages ~5- 
10 Myr in order to avoid the complications introduced 
by extinction, ongoing accretion, and veiling. The age 
sequence of ~M6 objects necessarily contains younger 
objects that show signs of accretion activity and higher 
extinction. The majority of our targets have Ay~0 with 
no evidence for ongoing accretion, although a few objects 
may retain cool disks. The majority of objects mak- 
ing up the 5-10 Myr sequence are confirmed low-mass 
members of the Upper Scorpius OB association, with 
additional young brown dwar fs from the ~10 Myr-old 
TW Hydrae Association (e .g. iZuckerman fc Son^ 120041 : 
IBarrado v Nava scues''2006y The sample properties are 
summarized in Table [21 

2.2. Observations 
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Observations were made between 2000 December and 
2009 April as part of the NIRSPEC Brown Dwar f 
Spectroscopic Survey fBDSS. lMcLean et al.|[2003ll2007f) . 
NIRSPEC is the cryogenic cross-dispersed echelle spec- 
trometer on_thfi_JCeckIIJ^ni__te^ on Mauna Kea, 
Hawaii ( [McLean et al. l fl99a [2000l) . NIRSPEC has two 
modes, a cross-dispersed echelle mode and a non-echelle 
mode in which the echelle grating is replaced with a 
mirror and the cross-disperser provides the spectral dis- 
persion. Both of these modes were employed with the 
NIRSPEC-3 (N3) filter, which approximates standard J- 
band coverage (1.143-1.375 /im). In echelle mode eight 
usable dispersion orders (65 to 58) are captured on the 
detector. Because the spectral interval captured by the 
detector is slightly smaller than the free spectral range 
in each order, there are small gaps, increasing with wave- 
length, in the total spectral coverage. The exact wave- 
length ranges for each dispersion order are listed in the 
headings in Setion 14.11 The slit width was two pix- 
els (0'.'38) for non-echelle observations and three pixels 
(G'.'432) for echelle observations. The resolving power 
in J-band is approximately R=A/AA=2,000 ("medium" 
resolution) in non-echelle mode and R=20,000 ("high" 
resolution) in echelle mode. Throughout the paper 
medium-resolution J-band spectra will be refered to by 
the filter name "N3" and high-resolution spectra by the 
number of the dispersion order. 

Observing methods are described in detail in 
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120031 (medium resolution) and 
120071 (high resolution); the following 
is a brief summary and explanation of departures from 
those methods. In both modes observations were made 
in pairs, nodding along the slit between each observation 
so that traces were separated by 20" on the 42"-long 
slit (medium resolution) and by 7" on the 12"-long 
slit (high resolution). For some recent high-resolution 
observations the size of the nod was reduced to avoid 
problems from an intermittent quadrant in the slit- 
viewing camera. During these occasions (2006 May 
and later) the nod size was at least 2" so that the 
dispersed traces would be well-separated on the slit. 
Integration time was 600 seconds per nod for all echelle 
mode observations except those on 2006 May 19 and 
20, which were 300 seconds per nod. Medium-resolution 
observations were typically 300 seconds per nod. Total 
integration times per object are listed in Table [TJ AO V 
stars were observed at an airmass very close to that 
of the target object (typically <0.1 airmass difference) 
to calibrate for telluric absorption features. Arc lamp 
spectra were obtained at least once per night, and 
white-light spectra and corresponding dark frames were 
obtained for flat-fielding. 

2.3. Data Reduction Methods 

All of the observed data were reduced with the 
REDSPEC ID L-based software package'', described in 
[McLean et all (2003, 2007). The package performs stan- 
dard bad pixel interpolation, dark subtraction, and fiat- 
fielding as well as spatial rectification of curved spectra. 
Spectra are rectified and extracted in differenced nod 
pairs so that the sky background and OH emission lines 
are removed. Spectra were extracted by summing over 

^ http:/ /www2. keck. hawaii.edu/inst/nirspec/redspec. html 
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7-15 pixels dependent on seeing then subtracted again 
to produce a positive spectrum with residual sky emis- 
sion features removed. AO V star spectra were reduced 
the same way, and target spectra were divided by the 
AO V star spectrum to remove telluric features and mul- 
tiplied with a Teff — 9500 K blackbody spectrum to 
restore the true slope. At medium resolution Pa/3 ab- 
sorption at A= 1.2822 /im in the AO V star spectrum is 
removed by interpolation, and at high resolution order 
59 was not divided by an AO V as it lacks strong tel- 
luric lines. Medium-resolution spectra were wavelength 
calibrated using Ne and Ar arc lamp lines, and high- 
resolution spectra with OH night sky lines, which were 
found to be highly stable and well-distributed across or- 
ders. Seven high- resolution dispersion orders (58-65, ex- 
cepting 60) were reduced. Order 60 was not reduced be- 
cause the OH night sky lines are blende d with O2 emis- 
sion bands at 1.26-1.28 /im (.Rousselot et al.ll2000f) mak- 
ing wavelength calibration and sky subtraction consider- 
ably more difficult. Each reduced spectrum was contin- 
uum normalized, and multiple nod pairs were averaged 
together to increase SNR. High-resolution spectra were 
shifted to the heliocentric reference frame. 

Because the sample objects were selected in order to 
minimize the effects of extinction, only GY 5 required 
correction for its relatively large extinction, Ay~5. The 
medium-resol ution N3 spect rum was dereddened as de- 
scribed in M cGovernl (120051 ). High-resolution observa- 
tions were not dereddened because the spectra cover such 
narrow wavelength ranges and are considered individu- 
ally in the analysis, making differences in the pseudo- 
continuum level caused by reddening relatively unimpor- 
tant. 

3. MODEL ATMOSPHERES AND SYNTHETIC SPECTRA 
3.1. The PHOENIX Code 

The PHOENIX code was developed by Peter Hauschildt 
and collaborators beginning about 20 years ago, origi- 
nally to model radiative transfer in supernova remnants. 
More recent versions of the code have been successfully 
applied to m odeling stellar, subste l lar, and planetary 
atmospheres jH auschil dt et all 119971 11999 '. AUar d et all 
120011: iBarman et al. 200W^ Some of the capabilities an d 
features of PHOENIX are reviewed bv lBaron et al.l (|2003f ). 
Atmospheres are constructed by calculating tempera- 
ture and pressure in each of 64 spherically-symmetric 
layers. The PHOENIX code determines radiative fluxes 
through the layers under the assumptions of hydrostatic 
equilibrium, radiative-convective equilibrium, and chem- 
ical equilibrium. The code begins with the results of 
a previously calculated structure, as an initial guess, 
and iterates toward the new solution by changing the 
radiative-convective temperatures at each layer via a 
modified Unsold- Lucy t emperature correction procedure 
(jHauschildt et al]l2003f l. Models were run for at least 
20 iterations, which was sufficient in most cases to reach 
convergence. Convergence of a model is based on energy 
conservation and determined by an upper-limit on the 
percentage difference between the radiative-convective 
flux and the desired model flux {a T^tfj), typically 5% 
or less at all layers. The output for converged atmo- 
spheres includes spectral energy distributions (SEDs) 
from 10 A to 1 mm with a resolution of 4 A through the 



optical and near-infrared. From a converged atmospheric 
structure the PHOENIX code is used to produce a syn- 
thetic spectrum for a selected wavelength range (i.e., the 
J band) and numerical resolution, up to 500,000 wave- 
length elements in total. 

3.2. Brown Dwarf Model Atmospheres 

To match our observational sample we calculate model 
atmospheres for effective temperatures from 1800 K to 
3000 K in intervals of 50 K and surface gravities from 
log(g) = 3.0 to 6.0 in intervals of 0.1 dex. For this tem- 
perature range the dusty version of the PHOENIX model 
code is applicable. In this version, dust is formed in 
chemical equilibrium conditions and remains where it 
is formed in the atmosphere; the removal of dust by 
dynamical processes is ignored. At low temperatures 
where more dust is created {Teff < 2200 K), the mod- 
els often do not converge because of the drastic ef- 
fects dust production has on the temperature structure. 
Instead, the temperature iterations will cause fluctua- 
tions around the converged solution. In these cases 
we manually chose the temperature-pressure structure 
with the smallest difference between calculated and pre- 
scribed flux (cr Teff), which was comparable to the dif- 
ferences for automatically converged models. All dust 
is created as pure, single-species, spherically-symmetric 
grain s that follow an interstellar grain-size distribu- 
tion (jAllard et al.l I2001I ). Model calculations incorpo- 
rate some updates to the [AUard ct al. (2001) models, 
including solar abund ances from lAsplund et al.l (|20051 ). 
the FeH l ine list fromlDulick et all (|2003f) . the CrH line 



list f rom [Burrows et al.l (|2002fl . and the H2O line list 



from lBarber et aLn200(J ). The incorporation of the new 
FeH line list causes the most drastic change from previ- 
ous generations of PHOENIX models and is discussed in 
detail in § 15.21 The differences are most evident in the 
medium- and high-resolution spectra and have little ef- 
fect on the atmospheric structures. We integrate under 
the SEDs and compare the total flux to a T'^jf , and the 
differences were typically less than 2% but up to 5% for 
models with Teff < 2000 K. Figure [T] shows the struc- 
tures and SEDs for a range of effective temperatures at 
a surface gravity of log(g)=5.0 and for a range of sur- 
face gravities at Te//=2600 K. The non-smoothness in 
the structures at cool temperatures and low pressures 
are induced by dust formation and are mainly above the 
near-infrared photosphere, which lies at approximately 
lO'^-lO^ dynes cm^^ depending on effective temperature 
and surface gravity. 

Pressure-broadening of strong lines impacts both the 
structure and the emergent spectrum of an atmosphere 
by adding opacity at wavelengths up t o thousands of 
angstroms outside of the line center (e.g. iBurrows et ahl 
120001: lAllard et all 1200 ID . The numerical treatment of 
Van der Wa als broadening i n the PHOENIX models is 
described in iSchweitzer et aTl (|1996D . A more detailed 
treatment of the line broadening incorpor ates multiple 
pertu rbers in calculating the line wings ( jAllard et al.l 
120051 ). We tested this treatment and found that the dif- 
ferences in structures, SEDs, and high-resolution spectra 
were minimal for the temperatures and surface gravities 
under consideration. We also experimented with how 
widely the K and Na line profiles are calculated from 
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the line center. This is particularly important for alkali 
lines in the optical that have wings extending thousands 
of angstroms and contribute a large amount of opacity; 
when the line profiles are truncated it affects the atmo- 
spheric structure. The minimum width of the line profile 
above which there was not a noticeable difference was 
used. The consideration of line broadening will be more 
crucial for cooler atmosphere models (i.e. T dwarfs) for 
which dust opacity no longer dominates. 

Medium-resolution synthetic spectra were calculated 
for the wavelength range 1.1-2.5 /im at 0.3 A resolu- 
tion to oversample the instrumental resolution for the 
non-echelle (medium resolution, R~2,000) mode of NIR- 
SPEC at all near-infrared bands. For the current anal- 
ysis only the J-band spectra are considered, but future 
papers will extend this study to H- and X-band BDSS 
data. High-resolution synthetic spectra were calculated 
at 0.03 A resolution for the wavelength range 1.1-1.4 /im 
to cover the entire J band and oversample the instru- 
mental dispersion of NIRS PEC by a f actor of 5-6 at 
those wavelengths (McLean 'etall 120071 Table 2). The 
output of the PHOENIX code can be used to identify the 
most prevalent atomic, molecular, and dust grain opacity 
sources. 

3.3. Spectral Fitting 

The MPFIT IDL code develeped bv lMarkwardtl (pOOOl ) 
was used to implement Levenberg-Marquardt least- 
squares minimization and determine the best-fit model 
parameters for each spectrum individually. The fit pa- 
rameters are T^ff and log(g) from the grid of atmosphere 
models. Projected rotational velocity {vsini) is incorpo- 
rated using standard Unsold-type profile and convolu- 
tion. Radial velocity (RV) is implemented as a shift in 
the observed spectrum onto the model wavelength array. 
A subroutine of the fitting code linearly interpolates be- 
tween calculated models for each flux value in a given 
wavelength range in order to draw from an essentially 
smooth model grid. Once a new model is created for T^f / 
and log(g), the usini kernel and RV correction are ap- 
plied. The usini kernel uses a limb-darkening coefhcient 
of 0.6, which is typical although not well constrained. 
The flux values of synthetic spectra created using differ- 
ent limb-darkening coefficients differ by less than a tenth 
of one percent on average even at the lowest surface grav- 
ities and highest wsini values considered. 

The MPFIT routine was implemented two different 
ways. First, the initial T^ff and log(g) parameters 
were randomly drawn from the continuous range cov- 
ered by linear interpolation of the calculated synthetic 
spectra. Some of the resulting best-fit parameters had 
much higher values than was typical, indicating the 
code had found a local minimum rather than a global 
minimum, which is a known drawback of Levenberg- 
Marquardt least-squares minimization. The best-fit re- 
sults were not sensitive to initial values of vsmi and RV. 
The best-fit parameters with with lowest values were 
averaged and set as the initial parameters for the second 
implementation of MPFIT. Following the Monte Carlo 
method, the initial guesses were kept constant and each 
flux point of the observed spectrum was randomly resam- 
pled from within the noise on that pixel. The distribu- 
tion of the best-fit parameters then indicates the uncer- 



tainty on the results from the noise in the data, which are 
adopted as our formal (relative) uncertainties, knowing 
that uncertainties in the model are much more difficult to 
charact erize. Thus, like similar studies bv lMohantv et aTl 
(j2004aD . for example, our uncertainties represent rela- 
tive precision and not accuracy. A veracious description 
of the uncertainties in the atmosphere models would re- 
quire, for example, comparisons of res ults from different 
sets of models (e.g. iHelling et al.ll2008[ ). which is beyond 
the scope of this paper. Typical results of this procedure 
are illustrated in Figure [D All spectra were fit using 
10** iterations of each implementation of MPFIT. Sev- 
eral spectra were fit with both 10'^ and 10** iterations, 
and the differences in resulting best-fit parameters were 
negligible. 

The results for each spectral fit are presented in Ta- 
bles [SJ^fTl Tables [3] and U] present the best-fit parame- 
ters for effective temperature and surface gravity, respec- 
tively, of each N3 spectrum and high-resolution disper- 
sion order for the spectral type sequences. Each result 
is coded according to the quality for all orders based on 
visual inspection. No annotation corresponds to a good 
fit with minimal mismatches between the observed and 
synthetic spectra. Poorer fits are denoted by (value) for a 
few mismatches in depth or wing-shape for the strongest 
lines and by -{value)- for several mismatches. Entries 
that are struck are very poor fits that were not used in 
determining the adopted best-fit values. Note that the 
annotations arc the same for both tables, even though it 
might be the case that the gravity is acceptable and the 
temperature is causing the poor fit, or vice versa. The 
quality of fits for each order are described in 14. II and for 
each parameter in 14.21 Tables [5] and [6] present the best- 
fit parameters for the ^M6 objects and are annotated 
the same way. The adopted Teff and log(g) values are 
the mean of results from each dispersion order, weighted 
by the sum of the squared errors for all orders. With 
this method parameters that result in a good fit for one 
order but a poorer fit for other orders receive less weight. 
The velocity results presented in Table [7] are determined 
by the same weighted-mean method. The results are dis- 
cussed by parameter in § 14.2. llH?^?^ The N3-fit results 
were not used in the determination of the adopted best- 
fit values. 

4. ANALYSIS 

In the following sections we discuss the medium- 
resolution (R'^2,000, referred to as N3) and high- 
resolution (R~20,000, referred to by dispersion or- 
der) synthetic and observed spectra in terms of major 
temperature- and gravity-dependent features. Then we 
summarize the fitting results for each parameter {T^ff^ 
log(g), vsini, and RV). 

4.1. Overview of Spectral Fits 

N3 (1.143 - 1.375 fim) — The N3 filter on NIRSPEC 
corresponds to the J band, the wavelength regime in 
which flux of M dwarfs peaks. We use "N3" as short- 
hand for medium-resolution J-band spectra throughout 
the following sections. The morphology of observed spec- 
tra for M dwarfs with different spectral types and ages is 
shown in Figure[31 The overall shape of the N3 spectrum 
is sensitive to effective temperature, while the strengths 
of atomic lines are sensitive to surface gravity. At lower 
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temperatures the pseudo-continuum shape becomes more 
dependent on surface gravity. For most objects the N3 
fits are very good, with the exception of the region be- 
tween 1.20 and 1.24 /xm, which is Ukely missing FeH 
opacity in the atmosphere models (see § I5.2p . and at 
1.28 /im where young and late-type objects might have 
weak Pa/3 emission. For some objects the depth of the 
H2O band starting at 1.335 /im is not well matched, with 
the band being too strong in the atmosphere models for 
later- type objects. The apparent mismatch in the H2O 
depth could also be a results of a poor fit to the pseudo- 
continuum level just before the H2O band. Qualita- 
tively the medium-resolution synthetic spectra provide 
good fits to the observed spectra, but the best-fit pa- 
rameters are often discrepa nt from values predict ed by 
the evolutionary ni odels of iChabrier et al.l ()2000t ) and 
iBaraffe et afl pOOl . 

Order 65 (1.16496 - 1.18207 fim) — The shortest- 
wavelength dispersion order features a strong K i dou- 
blet (resolved into a triplet at high resolution), a few 
weak Fe and Ti lines, and molecular lines of FeH and 
H2O. The K I lines are sensitive to both temperature 
and surface gravity because both decreasing temperature 
and increasing surface gravity contribute to pressure- 
broadening (Figure [3]) . The strength of the molecular 
lines can break the degeneracy between high gravity and 
low temperature, particularly the strongest H2O line at 
1.16579 /im, which is weakly sensitive to gravity but 
strengthens greatly with decreasing temperature. The 
synthetic spectra reproduce the pressure-broadening of 
the K I lines for later spectral types (decreasing effec- 
tive temperatures) and increasing age (increasing surface 
gravity) remarkably well. One difference between the or- 
der 65 observed and synthetic spectra is the shape of the 
pseudo-continuum between the K i lines at low tempera- 
tures. In the synthetic spectra the pseudo-continuum is 
more rounded as opposed to slightly flattened in the ob- 
served spectra. This is likely an indication that the alkali 
line profiles are not calculated as far into the wings as 
they should be. This was a known problem in the optical 
where the alkali lines are so strong that the calculation 
of the line profiles affects the atmospheric structure, but 
it is only a minor effect in the near-infrared. 

The atmosphere models provide excellent fits to the or- 
der 65 observed spectra because of both the well-modeled 
K I lines and the lack of strong FeH lines in this particular 
wavelength region. Although even the strongest H2O line 
is not strong enough to break the degeneracy between 
temperature and gravity with our fitting method, the 
linear relationship is so tight (see Figure [H first panel) 
that constraining either temperature or gravity based on 
spectral type or evolutionary models, respectively, allows 
only a narrow range of values for the other parameter. 
In general the best-fit parameters from order 65 provide 
excellent fits to the other orders (with the notable ex- 
ceptions of 62 and 63, see below), but in some cases the 
model fits produce K i lines that are too shallow in or- 
der 61 and Al i lines that are too deep with wings that 
are too broad in order 59. This may however be a conse- 
quence of the underrepresented FeH opacity in the atmo- 
sphere models that is changing the relative level of the 
pseudo-continuum and not a failure of the best-fit model 
parameters. 

Order 64 (1.18293 - 1.20011 jim) — Order 64 con- 



tains several atomic lines of Fe and Ti that are insensitive 
to temperature and gravity except at the low tempera- 
ture/high gravity region of the parameter space covered 
in this analysis. This order provides less consistent fits 
for all but the earliest spectral type young objects be- 
cause the FeH lines from the 0-1 band of F^ A-X"' A 
at 1.1939 /im that strengthens at lower temperatures 
and higher gravities are not reproduced by the atmo- 
sphere models (see § I5.2[) . Weaker atomic lines of Mg 
(1.1834 /im) and Ti (1.1896 & 1.1953 /xm) identified in 
the model spectra are present in the observed spectra, 
although they are more blended with other features in 
the observed spectra than the atmosphere models predict 
(Figure [5]). Order 64 produces good fits for the earliest 
spectral type young objects, with the fits becoming worse 
with later spectral types, resulting in high temperatures 
and low surface gravities that produce too- weak K i lines 
in order 65 and 61 and too-strong Al i lines in order 58. 

Order 63 (1.20168 - 1.21938 fim) — The only atomic 
lines in order 63 are very weak so that the spectral fea- 
tures are primarily FeH and H2O. As mentioned above, 
the FeH lines in the atmosphere models are too weak 
for all but the earliest spectral types of young objects. 
Figure [B] shows the general correspondence between ab- 
sorption features and the increasing mismatch in their 
strengths at older ages and later spectral types. In both 
the models and the observations the strength of the FeH 
lines increases more dramatically with decreasing tem- 
perature and later spectral type than with decreasing 
gravity (the latter is more apparent at even lower tem- 
peratures than are shown in Figure^. For young objects 
the correspondence between the locations of the lines, if 
not the strengths, is enough to provide a reasonable ra- 
dial velocity measurement in most cases (see § I4.2.3P . 

Order 62 (1.22093 - 1.23899 fxm) — Order 62 is also 
dominated by FeH lines that are not strong enough in the 
atmosphere models (Figure [7]), particularly strong lines 
from the Q-branch of the F^A7/2-X''A7/2 system from 
1.221 to 1.223 Jim and P- and j?-b ranch lines through- 
out the order (jMcLean et al.l [20071 ). Therefore order 62 
is very similar to order 63 in terms of fitting observed 
spectra with synthetic spectra: for the most part the fit 
parameters are inconsistent except for radial velocity. 

Order 61 (1.24081 - 1.25913 fim) — The longer- 
wavelength K I doublet in the J band falls in order 61 
(Figure [5]) . While the atmosphere models again repli- 
cate the dependence of the K i lines on temperature and 
gravity very well, the weaker lines in this order are more 
problematic than in order 65 for two reasons. First, the 
model predicts evenly-spaced TiO lines that increase in 
strength with decreasing temperature, particularly at low 
gravity, but such features are not seen in the observa- 
tions. This is likely because TiO lines in the observed 
spectra are blended with FeH lines that are too weak 
in the atmosphere models. Second, the lack of strong 
FeH in the models causes the wings of the K i lines and 
the pseudo-continuum shape to not match the observa- 
tions. Therefore order 61 provides less consistent model 
fits than order 65 despite the strong, well-modeled K i 
lines. 

Furthermore, the best-fit parameters do not fit other 
orders as well as the fits from order 65 that are also 
anchored by K i. For all but the earliest spectral types 
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(M6-M7), the best-fit parameters from order 61 result in 
Al I lines that are too strong in order 58 for the young ob- 
jects and for the field objects, an unrealistic combination 
of high temperature and high gravity or low temperature 
and low gravity, both of which arc manifested as poor fits 
in the other orders. 

Order 59 (1.28262 - LSOlSl^xm) — The atomic lines 
in order 59 are weak, but they have interesting depen- 
dence on temperature and gravity. The strongest lines 
are Ti and Mn, with Fe and Cr lines being too weak 
to distinguish from noise on the observational spectra. 
The strongest Ti lines (1.2835 and 1.2850 /im) and Mn 
line (1.2903 /xm) become broader and shallower with de- 
creasing temperature at high gravity in the atmosphere 
models, but the trend is not as apparent in the observed 
spectra of field objects (Figure [5]). At low gravity the 
lines do not change much with temperature. At high 
gravity and low temperature order 59 also contains FcH 
lines that are not strong enough in the models, render- 
ing this order poor for model fits. The narrow lines from 
1.297 to 1.301 /im in the spectra of the young objects are 
weak telluric lines that were not removed because this or- 
der contains Pa/? in the spectrum of an AO V star. Order 
59 provides consistent fits only for the earliest spectral 
type young objects. 

Order 58 (1.3044? " 1.32370nm) — The longest- 
wavelength order contains an Al doublet that, in the 
atmosphere models, is sensitive to temperature at high 
gravities but nearly unchanged with temperature at 
log(g) < 4.0 (Figure [TU]). The wings of the fines 
cfiange shape with decreasing temperature as they be- 
come blended with strengthening molecular lines. Tfiis 
behavior is also seen in the high-gravity field objects, but 
not the young objects, for which the Al lines are stronger 
in the atmosphere models than in the observations. This 
is likely a consequence of the ubiquitous FeH absorption 
lowering the pseudo-continuum level making the Al lines 
appear relatively weaker than in the models, in which the 
FeH is too weak. Order 58 provides only slightly more 
consistent fits than order 59. 

4.2. Spectral Fitting Results 
4.2.1. Effective Temperature 

Tables [3] and [5] present previously determined effective 
temperatures from the literature (column 3), the effec- 
tive temperatures derived from model fits to medium- 
resolution N3 spectra (column 4) and each of seven 
dispersion orders of high-resolution spectra (columns 
5-11), and our adopted best-fit temperatures (column 
12). Previously published effective temperatures are a 
combination of measurements of bolomctric luminosity 
(e.g. IVrba et all 120041 : IGolimowski et al. 2004 ), spectral 
type to temperature c onversions ( Luhman et al.l 120031 : 
iMohantv fc Basrill2003D . and com parison with syn thetic 
spectra at medium resolution fe.g. lMcGovernI [20051 ). For 
many of the objects the range of temperatures reported 
in the literature is 300-400 K. Almost all of our adopted 
temperatures fall within or are very close to the range 
of previously determined temperatures; notable excep- 
tions are TWA 5B and the field objects, discussed below. 
Adopted temperatures, which are from high-resolution 
fits, are also very similar to best-fit results from the N3 
spectra, within ~60 K for objects earlier than M7 with 



the exception of Gl 406, also discussed below. The re- 
sults become less consistent for the later spectral type 
objects, with the spectral fits producing higher tempera- 
tures than predicted by spectral-type effective tempera- 
ture relationships and bolometric measurements. For the 
most part the adopted temperatures decrease with later 
spectral types and are within ~100 K for objects of the 
same spectral type, again with the notable exception of 
TWA 5B. Both the N3 and echelle spectra of TWA 5B 
are lower signal-to-noise than average for our observa- 
tions and the spectra are likely contaminated by light 
from the J=7.67 mag primary ~2.5" away. Adopted 
effective temperatures for the field dwarf sequence also 
decrease for later spectral types, but they are even more 
dissimilar to previous temperature measurements. 

The measured effective temperatures are most dis- 
crepant for the latest spectral type young objects, in par- 
ticular TWA 5B, and the field (old, high gravity) objects 
Gl 406, LP 402-58, LP 412-31, and 2MASS 0140-^27. 
For all of these objects the temperatures measured via 
spectral fitting are 200-350 K hotter than the effective 
temperatures determined by spectral type-effective tem- 
perature scales and bolometric luminosity measurements. 
The higher temperatures measured from spectral fitting 
are likely a consequence of the increasing importance of 
FeH absorption for cooler and higher gravity objects that 
is not reproduced by the atmosphere models (see i; 15. 2p . 

The age sequence of ~M6 objects (Table [S]) shows a 
very small spread in measured effective temperature. In- 
terestingly the hottest temperature was measured for the 
unresolved binary Gl 577BC, which is the earliest spec- 
tral type object in the subsample. The average tempera- 
ture for the objects is 2880 K, with a standard deviation 
of 60 K, which, although cooler than the temperature 
of 2990 K for a young M6 on the Luhman et al. (200$ 
scale, is promising for the development of a spectral type- 
effective temperature scale that is consistent for objects 
with a range of properties (e.g., age, rotation, and bi- 
narity) and calibrated by benchmark objects for which 
observed properties are successfully reproduced by atmo- 
sphere models. 

Uncertainty resulting from noise in the observed spec- 
tra is determined via the standard deviation of the best- 
fit effective temperatures from flux-re-sampled spectra, 
which is typically < 10 K but as high as 30 K. Systematic 
uncertainty in the accuracy of the atmosphere models 
likely dominates the total uncertainty in the results, al- 
though as discussed above uncertainties in the models are 
difficult to quantify (see § I3.3P . The standard deviations 
of best-fit effective temperatures from 10"* iterations with 
varied initial parameters are typically 100-200 K. This is 
similar to the standard deviation of the best-ft effective 
temperatures from all orders, which is additionally sub- 
ject to varying levels of fidelity between the observations 
and the models (described in ij |4.ip . The systematic un- 
certainties in effective temperature are much lower for 
the medium-resolution spectra, typically < 50 K. The 
implications of the relative uncertainties from high- and 
medium-resolution spectra are discussed in § 15.11 

4.2.2. Surface Gravity 

The surface gravity measurements are presented in 
Tables |4] and [6] Previous studies have discovered 
gravity-sensitive features and spectral indices in the 
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near-infrared (e.g . iGorlova et al.1 120031 : iMcGovern et alj 
[2OOI lAllers et al.ll2007f) with low (R-300) and medium 
(R~2,000) spectra, notably a sharply peaked if-band 
spectrum (from increased H2O absorption shortward of 
the peak in H-hand flux) and weaker atomic lines (K i 
and Na i) . Our analysis differs from previous work in that 
surface gravity and effective temperature are determined 
simultaneously and at higher resolution. For young 
objects, the adopted log(g) from high-resolution fits is 
within 0.3 dex of the best-fit results for the medium- 
resolution spectra, with the exception of the TW Hydrae 
members and 2MASS 0608-27. For these objects the 
medium-resolution fits produced surface gravities nearly 
an order of magnitude lower than expected for their tem- 
perat ures and ages based o n the evolutionary models 
of C habrier et al. (2000) and lBaraffeeFaLl (pOOl ) (here- 
after DUSTYOO), an indication that high- resolution spec- 
tra are better for examining gravity-sensitive features 
than medium-resolution spectra. For field objects the 
medium-resolution fits produce gravities that are higher 
than physically reasonable according to the DUSTYOO 
evolutionary models. Such high gravities are likely a re- 
sult of the N3 fits being affected by the missing FeH 
opacity in the atmosphere models. The surface gravity 
measurements for Upper Scorpius objects have a mean 
of 3.87, which is very close to the mean of log(g)=3.91 
for 5 Myr objects between 2500 and 3000 K predicted by 
the DUSTYOO evolutionary models. As shown in Fig- 
ure 1181 the distribution of measured surface gravities is 
only slightly wider than predicted by the DUSTYOO evo- 
lutionary models, with outlying values for USco 66 AB 
(log(g)=4.26) and DENIS 1619-24 (log(g)=3.49). The 
relatively high wsim of USco 66AB or higher-order unre- 
solved binarity may contribute to the high surface grav- 
ity Fits to spectra of DENIS 1619-24, on the other 
hand, result in a low surface gravity despite the high 
usini and possible binarity. There is no straightforward 
explanation for the outliers, but they are consistent with 
the DUSTYOO evolutionary models within the systematic 
uncertainty in log(g), ~0.5 dex (see below). 

The surface gravity of field dwarfs as measured by fits 
to the high resolution spectra are somewhat higher than 
allowed by the DUSTYOO evolutionary tracks, but the 
measurements also have higher uncertainty because of 
the increasing mismatch between observed and model 
spectra for later spectral types and cooler temperatures, 
caused in large part by the atmosphere models lack- 
ing FeH opacity. The surface gravities measured from 
the N3 spectra are even higher, again suggesting that 
high-resolution spectra are required to accurately mea- 
sure surface gravity. It should be noted, though, that 
our adopted gravities are the weighted mean of results 
from individual orders that are often either much lower 
or much higher than what is predicted by the DUSTYOO 
evolutionary models. 

The age sequence of M6 objects (Table [5]) generally 
shows the predicted trend of increasing surface gravity 
with age, although it is not strictly monotonic. The sur- 
face gravities are within ~0.2 dex of what is predicted by 
the isochrones for their age and effective temperature. A 
notable exception is 2MASS 2234-)-40AB, which has the 
highest effective temperature and surface gravity of the 
young (<100 Myr) objects. However, the adopted tem- 
perature and gravity are driven upwards by the poor fits 



from orders 62 and 63, and temperature and gravity are 
degenerate such that lowering them both still produces 
a good spectral fit at high resolution. 

Uncertainty resulting from noise in the observed spec- 
tra is determined via the standard deviation of the best- 
fit surface gravity from flux-re-sampled spectra, which is 
typically < 0.1 dex but as high as 0.16 dex. As with effec- 
tive temperature, systematic uncertainty in the accuracy 
of the atmosphere models likely dominates the errors in 
the results. The standard deviation of best-fit surface 
gravities from the varied initial parameters and the stan- 
dard deviation of the best-ft surface gravities from all 
orders were typically 0.4-0.6 dex. Systematic uncertain- 
ties in surface gravity from the medium-resolution spec- 
tra were on average ~0.4 dex. The implications of the 
relative uncertainties from high- and medium-resolution 
spectra are discussed in § 15.11 

4.2.3. Projected Rotational and Absolute Radial Velocities 

Results for radial velocity and projected rotational 
velocity (RV and wsini) are presented in Table [71 
We use the well-measured radial velocity of Gl 406 
(19±1 km s~^) to establish systematic uncertainties in 
our measurements. The radial velocity of GL 406 mea- 
sured via spectral fitting is 18.7 km s~^, which is con- 
sistent with the previously measured RV. Uncertainties 
are determined using the Monte Carlo implementation of 
MPFIT (see § 13. 3p by taking the weighted mean of stan- 
dard deviations of the re-sampled-flux results. Typical 
uncertainties are 1-2 km s~^. 

Radial and projected rotational velocities are measured 
for nine objects for the first time: five Upper Scorpius 
members, the Taurus object CFHT Tau 7, the young ob- 
jects 2MASS 0608-27 and Gl 577BC, and the field dwarf 
LP 402-58. The average RV of the eight Upper Scorpius 
members in our sample is —7 km with a dispersion of 
<2 km s^^, consistent with the values of —5 km s~^ and 
—6 km s~^ det e rmine d by iMuzerolle et al.l ()2003[ ) and 
[Kurosawa et al.l (|2006f ) for five and thirteen very-low- 
mass Upper Scorpius members, respectively. The mea- 
sured radial velocity of SCH 1622—19 is furthest from the 
average of Upper Scorpius members, but still consistent 
with cluster membership considering our estimated un- 
certainty of 1-2 km s~^ . The measured radial velocity for 
CFHT Tau 7 is the peak of the histogram of T aurus mem- 
bers compiled by iBertout fc Genoval (|2006f ). providing 
more evidence of Taurus membership for this recently- 
discovered object. RV results for TW Hydrae members 
are self-consistent with a mean of 8 km s""'^ and a stan- 
dard deviation of <2 km s~^; however, they are system- 
atically lower than previous results by 1-6 km s~^. The 
RV measurements from our analysis are co nsistent with 
other confirmed TW Hydrae members (e.g. [Tor re s et al.l 
[2003h . The measured RV = -6.2 km s"^ for Gl 577BC 
is very similar to the most recently measured RV of the 
primary, -6.5 km s"^ (|Nordstrom et al.ll2004f ). 

In order to test the results from the fitting procedure 
in which RV is one of four free parameters, we also mea- 
sure RV by cross-correlating observed spectra with syn- 
thetic spectra calculated using the other best-fit param- 
eters {Teff, log(g), and usini). The differences were typ- 
ically less than 1 km for individual orders, but up 
to several km s~^ for order 63, particularly for field and 
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late-type objects. The average RV indicated the cross- 
correlation peak for each order was within 1-2 km s~^ 
of the weighted mean RV from spectral fitting for all 
objects, indicating that precision is not sacrificed by in- 
cluding RV as a free parameter in the spectral fitting 
routine. 

The derived vsini values presented in Table [7| are ei- 
ther consistent with or higher than previously measured 
values. The high vsini values for slow rotators are an 
unavoidable limitation of the instrumental resolution, 
which can be estimated using the known slowly-rotating 
standard Gl 406 (usini < 2.9 km s^^). Our measured 
vsini of Gl 406 is 8 km s~^, indicating the lower limit 
of measurable vsini for the instrumental resolution. Two 
other objects have similar measured vsini values (DE- 
NIS 1605-24 and SCH 1612-20); therefore, these are 
considered upper limits on their wsini. For the remain- 
ing objects the measured vsini is similar to or slightly 
higher than vsini values reported in the literature with 
the largest discrepancy for the field object LP 412-31. 
Both LP 412-31 and 2MASS 0140+27 are poorly repro- 
duced by even the best-fit model parameters because of 
the lack of FeH opacity in the atmosphere models, which 
become a major source of spectral features for late M 
dwarfs. When the model spectrum does not match the 
observations well enough to constrain vsini, the fitting 
routine tends to return a high value of vsini because it 
flattens out the model spectrum to compensate for mis- 
matched features, thus minimizing the x^- However, it 
is unclear why the best-fit for on object would produce 
a much higher vsini when the previously published val- 
ues are at or below the instrumental resolution, as is 
the case for LP 412-31. LP 412-31 is a strongly flaring 
M dwarf with very st rong Ha emission and mag netic field 
for its spectral type (iR eid et al.ll2002tlStelzer et al . 20061: 
iSchmidt et"an [2007: Rc iners fc B asri 2007|); perhaps the 
high vsini measurement in the near-infrared is related to 
activity (the previously publish ed vsini was measured in 
the optical bv IReid et al.ll200^ . The estimated precision 
of the vsini measurements is ±5 km s~^ based on the 
typical standard deviation of measurements from each 
order. 

5. DISCUSSION 

5.1. Measuring the Physical Properties of Young Brown 

Dwarfs 

Model fits to medium-resolution spectra produce con- 
sistent estimates of effective temperature but only weakly 
constrain surface gravity. On the other hand, high- 
resolution fits can be misleading because of the degener- 
acy between effective temperature, surface gravity, and 
vsini, which all contribute to broadening of the strongest 
atomic lines. Weaker molecular lines can help break this 
degeneracy, but the imperfect correspondence between 
observations and atmosphere models for the molecular 
lines, particularly FeH, confuse this issue. Therefore, the 
physical properties of brown dwarfs are likely best esti- 
mated by a combining both medium- and high-resolution 
for simultaneous fitting. 

Representative observed and best-fit model spectra are 
presented in Figures [TTHni From a comparison of the 
best-fit spectra for the young (5 Myr) M6 object DE- 
NIS 1605-24 (Fig. [11]) and the field (> 1 Gyr) M6 dwarf 



Gl 406 fFig. [T^ . it is evident that both the atomic and 
molecular lines are stronger in the field object that in the 
young object. The best-fit synthetic spectra are a better 
match for the young object than for the field object, par- 
ticularly in the wings of the atomic lines. As described 
in § 13.21 this is possibly a result of the simplified line 
profiles that are currently implemented in the PHOENIX 
code. However, it is likely that shape of the line wings 
are also affected by blended molecular lines, especially 
FeH, that are not yet properly reproduced by the model 
(see 15. 2p . It is promising, however, that in the spectra 
of DENIS 1605—24 there is a substantial correspondence 
between weak molecular lines in the model, especially 
evident in orders 63-61. The same lines are present in 
the best-fit spectra for Gl 406, but the line depths are 
far too shallow in the model even for this relatively hot 
object. Note that the missing molecular opacity in the 
model is also apparent in the medium-resolution spec- 
trum of Gl 406 (last panel, gray) compared to the model 
spectrum with the weighted best-fit parameters from the 
high-resolution spectra. However, individual molecular 
lines can only be examined at high resolution. 

Figures [T3] and [T3] are representative of the fits for 
later spectral type objects at different ages (~10 Myr: 
2MASS 1139-31 and >1 Gyr: LP 402-58). For 
2MASS 1139—31 there are considerably more mis- 
matches in the best-fit spectra than there are for the 
slightly younger and earlier-type object DENIS 1605—24. 
While the strong K i lines and many of the weaker fea- 
tures in the order 65 and 61 of 2MASS 1139-31 are re- 
produced, the strongest hues in orders 64, 59, and 58 are 
too deep in the model spectra. This can be attributed 
to the missing FeH opacity that dominates the molec- 
ular lines at these wavelengths: the pseudo-continuum 
level in the synthetic spectra is too high because of the 
missing opacity, resulting in atomic lines that are too 
deep when the synthetic pseudo-continuum is lowered to 
the observed pseudo-continuum. In the order 62 and 63 
spectra the increasingly poor reproduction of the FeH 
bands at lower temperature and higher surface gravity is 
evident, although some features are coincident in wave- 
length. The best-fit parameters also reproduce the over- 
all continuum shape of the medium-resolution spectrum 
of 2MASS 1139-31 (last panel ofFig.[Tni), although there 
is missing opacity from 1.19-1.24 /im, similar to Gl 406 
in Fig. [T21 The sharp drop in flux at 1.24 fim is less 
pronounced when a more detailed calculation of atomic 
line profiles are implemented in the code fsee 13.21) . The 
spectral fit for field M7 dwarf LP 402-58 (Fig. also 
shows the success of the fit for the K i lines in order 65 
and 61 and the too-deep Fe i, Mg i, and Al lines in or- 
der 64, 59, and 58. Similarly, the FeH bands in orders 
63 and 62 are too weak in the model, likely contribut- 
ing to the mismatches in atomic line depths at shorter 
wavelengths as well. The pseudo-continuum slope of the 
medium-resolution spectrum is flatter in the observed 
spectrum than in the model, suggesting that a cooler 
effective temperature provides a better fit. An effective 
temperature 100-200 K lower flattens out the pseudo- 
continuum but also produces a too-deep H2O band at 
1.34 fim and atomic lines that are too broad at high res- 
olution, even in combination with lower surface gravity 
and vsini values. Improving the FeH opacity in the at- 
mosphere models is a likely first step to resolving this 
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issue. 

Comparing our adopted values of Teff and log(g) 
to isochrones of t he DUS TYOO evolu tionary models 
(jChabrier et al.l[200 0: Baraff e et al.ll2002f ) provides a con- 
sistency check and a means of evaluating our results. 
Even though the DUSTYOO evolutionary models are 
poorly calibrated for the youngest objects, they can rule 
out unphysical values. For example, objects 5 Myr old 
with temperatures from 1800 K to 3000 K all have log 
surface gravities between 3.81 and 3.96, objects older 
than 500 Myr in the same temperature range have log 
surface gravities between 5.15 and 5.40, and no object 
of any age or effective temperature has a surface gravity 
over log(g)=5.40. 

All of our adopted temperatures and surface gravities 
are consistent with predictions of the DUSTYOO evo- 
lutionary models. Evolutionary tracks can be used to 
infer the masses and ages of objects from their mea- 
sured temperatures and surface gravities. Adopted Teff 
and log(g) values for the 5-10 Myr subsample (including 
2MASS 0608-27) and additional objects with ages <1- 
3 Myr are compared to the DUSTYOO tracks in FigurefTSl 
Within the error bars all of the objects are consistent 
with the properties predicted for young (age < 20 Myr) 
brown dwarfs (mass < 72 Mj). Two slightly anoma- 
lous objects are DENIS 1619—24 (too low gravity) and 
USco 66 AB (too high gravity), which are discussed 
above. All of the objects are well above the deute rium- 
burmng mass limit (-13 M r-'Chab rier et al.l[200l . 

The relative systematic uncertainties on the best-fit ef- 
fective temperatures and surface gravities from high- and 
medium-resolution spectra indicate a promising method 
of inferring the physical properties of young brown 
dwarfs. Medium-resolution J-band spectra strongly con- 
strain effective temperature of mid-late M dwarfs, but 
not surface gravity. Neither can medium-resolution spec- 
tra be used to measure usini, which may influence surface 
gravity measurements. High-resolution spectra, particu- 
larly of the shorter-wavelength K i doublet in NIRSPEC 
dispersion order 65, strongly constrain surface gravity 
relative to effective temperature. Therefore, we antici- 
pate that using a measurement of effective temperature 
from medium-resolution spectra to constrain the effec- 
tive temperature for high-resolution spectral fitting will 
increase the accuracy of properties derived using high- 
resolution spectra. The measured range of wsim values 
does not correlate with age or known binarity. Improve- 
ments in the model atmospheres that more consistently 
reproduce the spectral features of both young and field 
brown dwarfs will improve the accuracy of vsini measure- 
ments and provide insight into the dependence of angular 
momentum evolution on mass and age. 

5.2. FeH Oscillator Strengths 

One of the glaring mismatches between observed and 
synthetic spectra is the apparently missing opacity from 
FeH bands, particularly in orders 62 and 63. Our model 
atmospher e calcul a tions incorporate new FeH line list 
from Du lick et all (I2003D as a n update to the previ- 
ously used iPhillips et af (|1987f l line list. However, one 
key difference between the use of these line lists in the 
PHOENIX cod e is that oscil l ator s trength scalings avail- 
able for the 'Phillips et al.' (1987) line list are not ap- 
phed to the .Duhck et al.. (,2003 ) line list. .Burgasser et all 



(I2003D show s that the absorption coefficients in the 
iDulick et al.l ([2003) line list correspond to absorption 
features in the medium-resolution J-band spectrum of a 
lat e-L dwarf (t heir Figure 2b). Our results indicate that 
the Dulick et al. (2003) line list, as incorporated into the 
PHOENIX code, fails to reproduce observed FeH absorp- 
tion features. 

Figures [15] and [T6| illustrate the effects of different 
FeH line lists and oscillator strength scalings on syn- 
thetic spectra. Figure [15] shows J-band spectra syn- 
thetic spectra (Te//=2600 K, log(g)=5.0, convolved to 
ma tch a medium-reso lution NIRSPEC spectrum) for 
the IDulick et al.l (|2003D FeH line l ist, the IPhillips et al.l 
dlMZD line hst, and the "scaled" IPhillips et al.l (|1987[ ) 
line list, for which red-optical and near-infrared FeH 
band strengths are multiplied by factors of 10-20 (and 
for one band, 150). It is evident from comparison with 
the observed spectrum of 2MASS 0140-H7 (M9) that 
the scaled FeH strengths using the IPhillips etall (119871 ) 
line lists match the observed spectra better, particularly 
the broad absorption features at 1.20 and 1.21 ^m that 
are essentially absent in spectra with unsealed line lists. 
At this resolution the difference between the unsealed 
line lists is imperceptible, but the high-resolution spec- 
tra provide further insight. In Figure [TBI the same syn- 
thetic spectra are convolved to match NIRSPEC disper- 
sion order 63 , and it is apparent that even the scaled 
IPhillips et al.l (|1987|) line list falls short of reproducing 
the observed spectrum of the M9 dwarf 2MASS 0140-^27. 
The Dulick et al.. (200 3) list provides more individual 
lines in this wavelength range, but the lines are still far 
too weak. In a future paper we will quantify the short- 
comin gs in the curr e nt FeH data by implementing scaling 
of the ' Duhck et all (|2003l ) line hst in the PHOENIX code. 

FeH is an important source of spectral features from 
the z band through the H band, par ticularly in the l onger 
wavelength half of the J band ([McLean et al.l [20031 : 
iCushing et an[2003l [2005t [McLean et al.ll2007f) . and the 
features strengthen with decreasing effective temperature 
until — L5 (Kirkpat rick et anil999( ). when FeH begins to 
weaken through the L-T transition then strengthen in 
mid-to-late T dwarfs, possibly as cloud-clearing allows 
flux to escape fro m lower in the atmosphere where FeH 
remai ns as a gas ([Burgasser et al.l [20021 : [McLean et "all 
2003; Cushing et al.l 120081 ). Therefore, further study of 
FeH in the atmosphere models is necessary to extend our 
near-infrared spectral fitting analysis to L and T dwarfs 
as well as to shorter and longer wavelength observations. 

Magnetic fields may also affect FeH strengths, and the 
Wing-Ford band of FeH (—1 /zm) has been used to mea- 
sure magnetic field strengths in M dwarfs. However, 
strong magnetic fi elds actually weaken th e magnetically- 
sensitive lines (see iReiners &: Basrill2007l Figures 2-10). 
Therefore, if magnetic fields need to be included in the 
atmosphere models in order to reproduce observed FeH 
bands, it is likely a second-order effect to the overall scal- 
ing of the oscillator strengths. 

5.3. Effects of Binarity 

Unresolved binaries in the sample are 
Gl 577BC, USco 66AB, 2MASS 2234-H40AB, and 
2MASS 1207-39Ab. These b inaries vary in mass 
ratio (g) fro m g ~1 (Gl 5 77BC: ILowrance et all [2005l. 
USco 66AB: IKraus et all [20051 2MASS 2234-h40AB: 
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Allers et all 120091) to q ~0.3 (2MASS 1207-39Ab: 
Mohantv et 31112007?) . Additional ly, DENIS 1619-24 i s 
a ca ndidate spectroscopic binary (jMohanty et al.|[2"005D . 
and iSlesnick et alJ (120061 ) suggests that SCH 1622-19 
might be an unresolved binary (a < 175 AU) based on 
overluminosity for its spectral type relative to other 
Upper Scorpius members on an H-R diagram. For our 
purposes the effects of unresolved binarity are minimal 
as long as the mass ratio is either very high or very 
low. The main concern is the effect of unresolved 
binarity on the estimate of vsini. Previous authors 
have noted that broadened absorption li nes may be 
indic ative of an unresolved binary (e.g. iSimon et al.l 
l2006f ). As shown in Figure [ITl the velocity shift between 
components of a close binary can mimic a higher value 
of wsim. However, the RV-shift needed to reproduce the 
spectrum of USco 66AB (wsim = 28 km s^^), using the 
DENIS 1605-24 (wsim < 7 km s"^) as a template, is 
at least several times larger than the maximum possible 
velocity difference of the known binary. Two possibly 
unresolved binaries, DENIS 1619-24 and SCH 1622-19, 
are among the highest vsini measurements, but they are 
not anomalous compared to other, apparently single, 
rapidly rotating objects. A future paper will explore this 
issue further by comparing spectral fitting results from 
spatially resolved (using NIRSPEC and Laser Guide 
Star Adaptive Optics, e.g. Konopacky et al. in prep.) 
and unresolved spectra of very low mass binaries. 

6. CONCLUSIONS 

We have estimated the effective temperature, surface 
gravity, absolute radial velocity and projected rotational 
velocity for a sample of young brown dwarfs and field 
M dwarfs. Comparison of high-resolution synthetic and 
observed spectra of young objects provide more accurate 
determination of surface gravity than from lower resolu- 
tion observations. However, the small wavelength range 
and lack of fidelity of FeH line strength in the atmosphere 
models result in anomalously high temperatures for late- 
M type objects. Despite the mismatch in line strengths, 
radial velocity is well-determined for all but the coolest 
field objects considered in this analysis. Additionally, the 
richness of lines at high resolution breaks the degeneracy 
between sources of line-broadening so that vs'mi is also 
well-determined for most objects. The principal results 
of the paper are as follows: 

(1) Adopted (weighted- mean) effective temperatures 
from high-resolution spectral fits for objects M7 and ear- 
lier fall within or very close to previously published esti- 
mates and results from medium-resolution spectral fits. 
For objects M8 and later, temperatures from medium- 
resolution fits are similar to previously published values 
and temperatures from high-resolution fits are systemat- 
ically too hot. 

(2) Adopted (weighted-mean) surface gravities from 
high-resolution fits neatly separate objects with Teff 
from 1800-3000 K according to age, with 1-10 Myr ob- 
jects lying broadly between log(g)=:3.5 and 4.3, while 
objects from 0.5-5 Gyr are in the log(g)=5.2-5.4 range. 
Surface gravities from medium-resolution fits are gen- 
erally similar to the adopted values except for M8-M9 
young objects, which are systematically low, and field 
objects, which are systematically high. 

(3) The age sequence of M6 objects yields an average 



temperature of 2880 ± 60 K. 

(4) For the Upper Scorpius objects the mean surface 
gravity is log(g)=3.87, which is very close to the value 
of log(g)=3.91 predicted by the DUSTYOO evolutionary 
model for 5 Myr objects between 2500 and 3000 K. 

(5) FeH features are not yet properly reproduced by 
atmosphere models, and this impacts fitti ng spectra to 
high-r esolution observations. The recent iDulick et al.l 
()2003D line lists provide better correspondence for indi- 
vidual features, but the lines are not strong enough in 
the models. Improving FeH in the atmosphere models is 
urgently needed. 

(6) Physical properties of brown dwarfs are likely best 
measured by comparing observed and synthetic spec- 
tra using a combination of medium- and high-resolution 
spectra simultaneously. 

This work represents the first in a series of analysis 
of medium- and high-resolution observed and synthetic 
spectra of brown dwarfs. Future work will extend the 
study to later spectral types, cooler atmosphere models 
(wh ere Tp f f and log(g) rn ight become more intertwined, 
e.g. IBurgasser et al.l l2OO60 . and different dust treatments 
(i.e. the PHOENIX- cond models). Further analysis will 
also explore effects of metallicity on atmospheric struc- 
ture and resultant spectra. While the focus thus far is 
on the large database of high-resolution spectra provided 
by the BDSS, some combination of high resolution and 
broad spectral coverage will likely prove ideal for accu- 
rately inferring the physical properties of brown dwarfs 
using atmosphere models. Therefore, we are developing 
a procedure to combine all available observed spectral 
data for simultaneous model fits. 
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Fig. 1. — Left: Temperature-pressure profile of model atmosphere structures calculated using the PHOENIX code. Filled circles mark the 
near-infrared photosphere (t)^—i 2 = §)• Right: Spectral energy distribution (SED) for the same model atmospheres. The top plots 
display structures and SEDs calculated at log(g)=5.0 for the entire range of effective temperatures. The bottom plots display structures 
and SEDs (offset by a constant) calculated at T^f f=2600 K for the entire range of surface gravities. Only half of the calculated structures 
and SEDs are shown for clarity. 
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Fig. 2. — The distribution of best-fit values found by the MPFIT procedure with varied input parameters (see §3.3) for the order 65 
spectrum (1.165 — 1.182 ^m) of DENIS 1619—24. Contours mark 10%, 50%, and 90% of the maximum number of points in 25 K, 0.05 
dex, 0.5 km ss^-'^ (RV) and 0.5 km ss^-*^ (tisini) bins. Filled circles mark our adopted best-fit values. Top left: Effective temperature and 
surface gravity show a strong degeneracy in this wavelength range. Top right: Radial velocity is well-determined independent of the other 
parameters. Bottom: Projected rotational velocity is slightly degenerate with both effective temperature and surface gravity. These plots 
do not show the value for each fit result, which is taken into account to determine our adopted best-fit parameters. 




Fig. 3. — Observed J-band spectra of young (age ~5-10 Myr) M dwarfs (left), field (age > 1 Gyr) M dwarfs (middle), and M6 objects 
with a range of ages (rigiit). Spectra arc normalized at 1.26 /im and offset by a constant. The left plot labels molecular absorption bands, 
the middle plot strong atomic lines, and the right plot the wavelength coverage of NIRSPEC dispersion orders. 
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Fic;. 4. — High-resolution (R~20,000) observed and synthetic spectra for NIRSPEC dispersion order 65. In both panels at top and 
bottom are observed spectra (SCH 1622-19: MS, 5 Myr; LP 412-31: M8, >1 Gyr; DENIS 1605-24; M6, 5 Myr; 2MASS 0608-27: M8.5, 
<100 Myr) whereas the five middle spectra are synthetic spectra with effective temperature and log surface gravity as labeled. Horizontal 
lines denote vertical offsets in the pseudo-continuum levels. The strongest absorption lines are marked by dot-dashed lines and labeled on 
the plot. The effects of increasing surface gravity (left) and decreasing temperature (right) on the pressure-sensitive K I lines are similar 
but can be disentangled at high resolution, e.g. using the strength of the H2O line at 1.16579 //m and other molecular lines. 
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Fig. 6. — The same as Figure|4]but for NIRSPEC dispersion orde r 63. The spectral features in this order are primarily H2O and FeH, 
the latter of which is too weak in the atmosphere models (see § I5.2| l . 
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Fig. 7. — The same as Figure|4]but for NIRSPEC dispers ion order 62. The spectral features in this order are primarily H2O and FeH, 
and again FeH is too weak in the atmosphere models (see S I5.2| l. notably the the Q-branch of the i<"*A7/2"^*^7/2 system from 1.221 to 
1.223 fim. 




Fig. 8. — The same as Figure l4l but for NIRSPEC dispersion order 61. Regularly-spaced lines in the synthetic spectra ar e fro m TiO. In 
the observed spectra TiO (if present) is blended with FeH, but FeH is considerably weaker in the atmosphere models (see ij 15.21 1. 
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Fig. 9. — The same as Figure |3] but for NIRSPEC dispersion order 59. 
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Fig. 10.— 



The same as Figure|3]but for NIRSPEC dispersion order 58. 
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Fig. 11.— Observed (gray) and best-fit synthetic spectra (black) for the young (5 Myr) M6 object DENIS 1605-24. The NIRSPEC 
dispersion orders are labeled, and the best-fit atmosphere model parameters are given in the bottom-right panel, which is the medium- 
resolution spectrum across the entire J band. 
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Fig. 12. — Observed (gray) and best-fit synthetic spectra (black) for the field (old, >1 Gyr) M6 dwarf Gl 406. The NIRSPEC dispersion 
orders are labeled, and the best-fit atmosphere model par ame ters are given in the bottom-right panel. Note the poor fit for the observed 
FeH features at both high- and medium-resolutions (see 8 15.2| l. 
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Fig. 13.— Observed (gray) and best-fit synthetic spectra for the young (~10 Myr) M8 TW Hydrae member 2MASS 1139-31. The 
NIRSPEC dispersion orders are l abele d, and the best-fit atmosphere model parameters are given in the bottom-right panel. Again, note 
the poor fit to FeH features (see §[32}. 



26 



Rice et al. 




1.170 



1.175 



1.180 



1.185 



1.190 



1.195 



1.200 



X 
3 



T3 

tsj 

E 




1.205 



1.210 



1.215 



1.225 



1.230 



1.235 





59 



1.245 



1 .250 1 .255 



1.285 1.290 1.295 1.300 



1.2 
1.0^ 
0.8 F 
0.6 
0.4 

0.2 
0.0 




58 



Teff = 2879 K 
log{g) = 5.31 
RV = -2.6 km s ' 
vsini = 17 km s"^ 



1.305 1.310 1.315 1.320 1.20 1.25 

Wavelength (|j,m) 



1.30 



Fig. 14. — Observed (gray) and best-fit synthetic spectra for the field (old, >1 Gyr) M7 dwarf LP 402-58. The NIRSPEC dispersion 
orders are labeled, and the best-fit atmosphere model parameters are given in the bottom-right panel. 
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Fig. 15. — Synthetic spectra produced using the same atmospheric structure (Tejj=2400 K and log(g)=5.0) but different FeH line lists 
and oscillator-strength scalings as indicated on the plot. The NIRSPEC medium-resolution spectrum of the M9 dwarf 2 MA SS 0140+17 
is plotted for compari son, and the horizo ntal bar marks the wavelength range of NIRSPEC order 63 shown in Figure 1161 While the 
difference between t he |Dulick et al.l {2003) and the unsealed Phillips e t al -i (1987) line lists are imperceptible at this resolution, the scaled 
IPhillips et al.l 1)19871 ') oscillator stre ngth s make a dramatic difference in the spectra from 1.205 to 1.240 fira, resulting in a much better 
match the the observed data. See l5.2lfor a detailed discussion. 
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Fig. 16. — Synthetic spectra produced using the same atmospheric structure (Tejj=2400 K and log(g)=5.0) but different FeH line lists 
and oscillator-strength scalings as indicated on the plot. NIRSPEC dispersion order 63 spectrum of the M9 dwarf 2MASS 0140+27 is plotted 
for comparison. The synthetic spectrum using the Dulick ct al.. (.2003 ) line list shows an improved correspondence over the i Phillips et al.l 
([ 1987 ) line list for individual features in the 2MASS 0140+27 spectrum, but the line strengths are too weak in the atmosphere model. See 
S l5.2l for a detailed discussion. 
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Fig. 17.— Observed NIRSPEC dispersion order 61 spectra of USco 66AB (top) and DENIS 1605-24 (bottom). The second-to-top 
spectrum are formed by shifting the spectrum of DENIS 1605—24 (log(g)=3.82, vshii < 7 km s^^) by 20 km s^^, the result added to the 
unshifted spectrum, and renormalizing. The second-from-bottom spectrum is the spectrum of DENIS 1605—24 convolved with a rotational 
velocity profile of 20 km s~^. Both the shifted and convolved spectra resemble the spectrum of USco 66AB (log(g)=4.26, vsim= 28 km s-i), 
indicating that binarity may influence the inferred -ysini and/or log(g). This is unlikely to be the case for USco 66AB, for which the known 
binary would induce no more than a few km shift in the spectrum (assuming a circular orbit), but could be a factor in as yet unresolved 
objects such as DENIS 1619-24 and SCH 1622-19. 
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Fig. 18. — Adopte d physical properties of young objects {log(g) < 4.3) plotted with the DUSTYOO tracks from lChabrier et al.l jgOOOl); 
IBaraffe et al.l II2002I ). Membership (and/or approximate age) are indicated by different symbols: circles for p Ophiuchi, Taurus, LHci233 
group, and a Orionis objects; squares for Upper Scorpius objects; triangles for TW Hydrae objects, and diamonds for candidate young field 
objects and companions. The error bars in the upper right represent typical uncertainties based on Monte Carlo simulations (see 13.31 1: 
systematic uncertainties in the atmosphere models arc likely larger. 
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TABLE 1 
Observing Log 



Object 


Spec. 


R.A. 


Dec. 




Int. Time 


UT Date 




Name" 


Type 


(J2000.0) 


(J2000.0) 


(mag) 


(seconds) 


of Observation'^ 


Ref. 


Gl 577BC 


M5.5 


15 05 49.9 


+64 02 50 


7.19 


1200 


2006 May 19 


1 












600 


2003 Mar 24 


3 


Gl 406 (Wolf 359) 


M6 


10 56 28.9 


+07 00 53 


7.09 


240 


2006 May 20 


1 










120 


2000 Dec 6 


2 


CFHT Tau-7 


M5.75 


04 32 17.2 


+24 22 11 


11.52 


1200 


2008 Mar 19 


1 












600 


2009 Nov 08 


1 


a Orionis 12 


M6 


05 37 57.5 


-02 38 44 


14.23 


6000 


2008 Dec 06 


1 












1200 


2003 Dec 4 


3 


UScoCTIO 66AB 


M6 


16 01 49.6 


-23 51 08 


12.91 


2400 


2007 May 29 


1 












600 


2003 May 12 


3 


DENIS-P J160514. 0-240652 


M6 


16 05 14.0 


-24 06 52 


12.84 


2400 


2007 May 29 


1 












1200 


2004 May 2 


3 


GY 5 (ISO-Oph 30) 


M6 


16 26 21.5 


-24 26 01 


12.57 


2400 


2007 May 31 


1 












600 


2003 May 14 


3 


2MASS J22344161+4041387AB 


M6 


22 34 41.6 


+40 41 39 


12.57 


2400 


2007 May 29 


1 












600 


2003 Dec 4 


3 


SCH J16121188-20472698 


M6.5 


16 12 11.9 


-20 47 27 


13.66 


4800 


2007 May 29 


1 












1200 


2009 Apr 7 


1 


LP 402-58 


M7 


23 36 43.9 


+21 53 39 


12.71 


3600 


2008 Dec 7 


1 












600 


2004 Nov 6 


3 


UScoCTIO 100 


M7 


16 02 04.3 


-20 50 43 


12.84 


2400 


2007 May 30 


1 












600 


2003 May 13 


3 


UScoCTIO 130 


M7.5 


15 59 43.7 


-20 14 40 


14.20 


4800 


2007 May 30 


1 












1200 


2004 Apr 30 


3 


LP 412-31 


M8 


03 20 59.7 


+18 54 23 


11.74 


2400 


2008 Dec 7 


1 












600 


2000 Dec 4 


2 


SCH J16235158-23172740 


M8 


16 23 51.7 


-23 17 26 


13.55 


3600 


2008 Mar 20 


1 












600 


2009 Apr 7 


1 


DENIS-P J161929. 9-244047 


M8 


16 19 29.9 


-24 40 47 


14.23 


2400 


2007 May 30 


1 












1200 


2004 Jul 22 


3 


SCH J16224384-19510575 


M8 


16 22 43.8 


-19 51 06 


12.35 


2400 


2008 Mar 19 


1 












600 


2009 Anr 7 


1 


2MASS J11395113-3159214 


M8 


11 39 51.1 


-31 59 21 


12.69 


3600 


2006 May 19 


1 












1200 


2001 Dec 30 


3 


2MASS J12073347-3932540AB 


M8 


12 07 33.4 


-39 32 54 


13.00 


2400 


2007 May 29 


1 












1200 


2001 Dec 30 


3 


2MASS J06085283-2753583 


M8.5 


06 08 52.8 


-27 53 58 


13.60 


4800 


2008 Dec 6 


1 












600 


2003 Dec 4 


3 


TWA 5B 


M8.5 


11 31 55.4 


-34 36 29 


7.67 


1200 


2006 May 20 


1 












600 


2002 Dec 24 


3 


2MASS J01400263+2701505 


M9 


01 40 02.6 


+27 01 50 


12.49 


2400 


2000 Dec 4 


4 












600 


2009 Nov 08 


1 


















References. — (1) this work. ('2')IMcLean 


et al.ll2003l. ("3") IMcGovernll2005l ('41IMcLean et al.ll2007l 





^ 2MASS, DENIS, and SCH object names are truncated in subsequent tables and in the text. 

^ From 2MASS All-Sky Point Source Catalog. The magnitudes for the close companions Gl 577BC and TWA 5B arc for 
the combined system because they arc unresolved in 2MASS. 

The first entry for each object refers to high-resolution (echelle) observations and the second entry, if present, refers to 
medium-resolution (non-echcllc) observations. References for these observations, if previously published, are given in the 
last column. 
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TABLE 2 
Sample Properties 



Object 


Spec. 






Name 


Type 


Notes 


Rcf. 


Upper Scorpius (~5 Myr) 


USco 66 AB 


M6 


q~l binary 


1, 2 


DENIS 1605-24 


M6 




3 


SCH 1612-20 


M6.5 




4 


USco 100 


M7 




1 


USco 130 


M7.5 


Estimated spectral type 


1 


SCH 1623-23 


M8 




4 


DENIS 1619-24 


M8 


Possible spectroscopic binary 


3, 4 


SCH 1622-19 


M8 


Possible unresolved binary 


5 


TW Hydrac (~10 Myr) 


2MASS 1139-31 


M8 




6 


2MASS 1207-39AB 


M8 


q~0.3 binary 


6, 7 


TWA 5B 


M8.5 


Companion 


8 


Other Young Objects 


GY 5 (ISO-Oph 30) 


M6 


p Ophiuchi, <1 Myr 


9 


2MASS 2234+40AB 


M6 


Lk Ha 233, ~1 Myr, q~l binary 


10 


CFHT Tau-7 


M5.75, M6.5 


Taurus, 1.5 Myr 


11, 12 


cr Ori 12 


M6 


(7 Orionis, ^3 Myr 


13 


Gl 577BC 


M5.5 


Companion, ~70 Myr, q~l binary 


14 


2MASS 0608-27 


M8.5 


Young Field, <100 Myr 


15 


I'iold Objects (>1 Gyr) 


Gl 406 (Wolf 359) 


M6 


RV standard 


16 


LP 402-58 


M7 




17 


LP 412-31 


M8 


Strongly flaring object 


17, 18 


2MASS 01404-27 


M9 




16 



References. — (l)IArdila et al]||200q). (2)IKraus et al.lll200H) . ( 3')IMartm et al.l l (200l) . (4) 
jMohantv et al.|[2005f) (5)ISlcsnick ct al. (200g). (6)'Gizis' (2002'). f7VMohantv ct al. f2007'). (8) 
ILowrance et al.Nl999D (9) Wilkiiie ct al. (1999), (10) AUcrs ct al. (2009), (li;Luliman I 200^ 
(12)IGuieu et al.ll|200q), (1 31lBeiar ct al., (|1999i') . fl 41 iLowrance ct al.i (l2005i'). f 151 iKirkpatricl^ 
i200g|) . (lel lMcLean et al.l I I2007I) . (17) IGizis et al.l I I2000I) . (IS HSchmidt et al.l I I2007I) . 



TABLE 3 

Spectral Fitting Results - Effective Temperatures 







Spec. 


Prev. 
















Adopted 


Object 




Type 


T f 


N3 


65 


64 


63 


62 


61 


59 


58 Te// 


USco 66AB 




M6 


2851-3000 


2973 


2927 


2939 


(2800) 


2250 


2988 


(3000) 


(2881) 2928 


DENIS 1605- 


24 


M6 


2850-3000 


2826 


2785 


2739 


(2537) 


(2800) 


(2977) 


2921 


2800 2794 


SCH 1612-20 




M6.5 


2630-2935 


2840 


2804 


2757 


2531 


2699 


(3000) 


2999 


2624 2777 


USco 100 




M7 


2672-2880 


2880 


2825 


2813 


-(2651)- 


-(2998)- 


(2912) 


(3000) 


(2746) 2849 


USco 130 




M7.5 


2550-2850 


2660 


2653 


2764 


(2800) 


(2883) 


(2826) 


(2978) 


(2741) 2805 


SCH 1623-23 




M8 


2455-2710 


2625 


(2489) 


(2780) 


(2351) 


-(2458)- 


(2741) 


-(2899)- 


-(2664)- 2636 


DENIS 1619- 


24 


M8 


2600-2710 


2726 


(2666) 


(2770) 




29Qg 


(2776) 


(2810) 


(2558) 2713 


SCH 1622-19 




M8 


2400-2710 


2594 


2594 


(2632) 


(2261) 


(2376) 


(2800) 


-(2939)- 


-(2792)- 2629 


2MASS 1139- 


-31 


M8 


2400-2710 


(2435) 


2486 


(2701) 


(2409) 


-(2900)- 


(2600) 


-(2837)- 


(2484) 2609 


2MASS 1207- 


-39AB 


M8 


2400-2710 


2485 


2602 


2694 


2405 


(2800) 


(2801) 


(2804) 


2486 2652 


2MASS 0608- 


-27 


M8.5 


2150-2555 


2200 


(2574) 


-(2585)- 


(2167) 


2276 


-(2900)- 


-(2727)- 


-(2474)- 2529 


TWA 5B 




M8.5 


2300-2550 


2360 


(2716) 


(2686) 


(2391) 


2787 


2g20 


(2904) 


(2814) 2731 


Gl 406 




M6 


2670-3058 


3000 


(3000) 


(2833) 




(2895) 


(2900) 


3QQQ 


ae©& 2911 


LP 402-58 




M7 


2600-2620 


2864 


(2950) 


-(2800)- 


(2900) 


(2799) 


(2900) 


3QQQ 


306& 2879 


LP 412-31 




M8 


2480-2638 


2919 


(2936) 


(2950) 


-(2273)- 


-(2805)- 


(2860) 


OgQQ 


OmO- 2806 


2MASS 01404 


-27 


M9 


2325-2400 


2824 


(2805) 


(2850) 


-(2215)- 


-(2846)- 


-(2900)- 


3QQQ 


374-4- 2753 


Note. — Sec 


§ 13.31 for cxplanat 


ion of annotations of the table entries. 













^ Effective temperature (or range of temperatures) determined by previous studies for the specific object or for an object of the same 
spectral type and similar age are from Basri et al. (2000); Dahn et al. (2002); Luhman et al. (2003); Gorlova ct al. (2003); Mohantv &z Basri] 
II2003I); [Golimowski et al.| (i200l); [Nloha nty ct al. (2004a); Vrba ct al. (2 004); McGove rn (.2005.) ; .Reiners Basri. (.2007.) ; .Slesnick ct al.„ 12006.') ; 
IScholz et aD (f2007h ; fslesnick et all JloOSl) ; ITeixeira et al.l j200gi l ; IWilking et al.l (fl99^ . 
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TABLE 4 

Spectral Fitting Results - Surface Gravities 



Object 


Agc^ 


Prev. 
Log(g)'= 


N3 


65 


64 


63 


62 


61 


59 


58 


Adopted 
log(g) 


USco 66AB 


5 Myr 


3.64-4.50 


4.21 


4.32 


4.20 


(4.69) 




4.72 


(3.80) 


(3.92) 


4.26 


DENIS 1605-24 


5 Myr 


4.25 


3.70 


3.60 


3.49 


(3.60) 


(4.47) 


(4.67) 


3.94 


3.14 


3.82 


SCH 1612-20 


5 Myr 


3.57-3.61 


3.50 


3.66 


3.46 


3.81 


3.53 


(4.67) 


3.80 


3.01 


3.68 


USco 100 


5 Myr 


3.55-4.50 


4.10 


4.09 


3.82 


-(4.70)- 


-(5.69)- 


(4.75) 


(3.72) 


(3.09) 


4.13 


USco 130 


5 Myr 


3.25-4.00 


3.82 


3.92 


3.56 


(4.70) 


(4.35) 


(4.67) 


(3.99) 


(3.12) 


4.05 


SCH 1623-23 


5 Myr 


3.27-3.32 


3.66 


(3.57) 


(3.48) 


(3.89) 


-(4.41)- 


(4.46) 


-(3.61)- 


-(3.00)- 


3.74 


DENIS 1619-24 


5 Myr 


4.00 


3.52 


(3.66) 


(3.38) 






(4.19) 


(3.21) 


(3.00) 


3.49 


SCH 1622-19 


5 Myr 


3.02-3.06 


3.70 


3.89 


(3.38) 


(3.79) 


(3.56) 


(4.77) 


-(3.85)- 


-(3.01)- 


3.76 


2MASS 1139-31 


10 Myr 


3.62-4.25 


(3.30) 


3.76 


(3.82) 


(4.03) 


-(5.59)- 


(4.30) 


-(3.57)- 


(3.33) 


3.98 


2MASS 1207-39AB 


10 Myr 


4.00-4.25 


3.30 


3.90 


3.50 


3.66 


(4.60) 


(4.65) 


(3.40) 


3.00 


3.83 


2MASS 0608-27 


<100 Myr 


4.00 


3.20 


(4.42) 


-(3.61)- 


(3.92) 


3.60 


-(5.57)- 


-(3.40)- 


-(3.00)- 


3.98 


TWA 5B 


10 Myr 


3.75-3.98 


3.10 


(4.12) 


(3.50) 


(4.20) 


4.34 




(3.60) 


(3.53) 


3.88 


Gl 406 
LP 402-58 
LP 412-31 
2MASS 0140-1-27 


>1 Gyr 
>1 Gyr 
>1 Gyr 
>1 Gyr 


5.15-5.40 
5.15-5.40 
5.15-5.40 
5.15-5.40 


5.67 
5.35 
6.00 
5.83 


(5.61) 
(5.55) 
(6.00) 
(5.51) 


(4.56) 
-(4.70)- 
(6.00) 
(5.10) 


3_g2 
(5.60) 
-(3.88)- 
-(3.61)- 


(4.85) 
(4.56) 
-(4.79)- 
-(4.88)- 


(5.60) 
(5.73) 
(6.00) 
-(6.00)- 


/[ QQ 


4^ 

/[ OQ 
3 /[Q 


5.19 
5.31 
5.53 
5.15 



Note. — See ^ |3.3| for explanation of annotations of the table entries. 
^ Age references: IGizis et aI1l2000| ; IPreibisch et al.|[200^ IZuckcrm an fc Song| |2004': 'Barrado y Navascuc^'2006'; 'Kirkpatrick"200^ 

Surface gravities determined by previous studies for young objects are from [Gorlova ct al. (2003); Mohanty ct al. (2004a); McGovernI (120051 ): 
IBarrado v Navascu^g ( 120061) ; 'Slcsnick ct al.' f200ff, '2008); Tcixcira ot al.' f2008). The range of surface gravity for field objects (age >1 Gyr) is 
set by evolutionary models of J3habnm^_ct_ (2000) and JBai\affc_ct ah (2002). 



TABLE 5 

M6 Sequence - Effective Temperatures 





Spec. 


Prev.^' 


















Adopted 


Object 


Type 




N3 


65 


64 


63 


62 


61 


59 


58 




GY 5 


Me'' 


2700-3050 


2871 


2939 


2814 


(2676) 


(2788) 


(3000) 


(3000) 


2845 


2876 


2MASS 2234-I-40AB 


M6 


2990 


2938 


2950 


2929 


(3000) 


(3000) 


3000 


2913 


2617 


2921 


CFHT Tau-7 


M6<^ 


2935-3024 


2740 


2805 


2801 


(2701) 


(2775) 


2946 


2969 


2760 


2825 


o-Ori 12 


M6 


2990 


2779 


2973 


2833 


2650 


2700 


3000 


2939 


2793 


2842 


DENIS 1605-24 


M6 


2850-3000 


2826 


2985 


2739 


(2537) 


(2800) 


(2977) 


2921 


2800 


2794 


Gl 577BC 


M5.5 


2900-292 


3000 


3000 


(2961) 


-(3000)- 


(2970) 


(2999) 


(3000) 


(3000) 


2988 


Gl 406 


M6 


2670-3058 


3000 


(3000) 


(2833) 


221^ 


(2895) 


(2900) 


3QQQ 


3QQQ 


2911 



Note. — See S 13.31 for explanation of annotations of the table entries. 
^ Effective temperature (or range of temperatures) determined by previous s tudies for the s pecifi c object or for an objec t 
of the s ame spectral type and similar age are from ILuhm an Rick? (T9W); "Luh man et~all ((2003); ILowrance eFall ( |2005|) : 
IMohant v et al. (2005); Wilking et al. ( 20051): IGatti et al.l (12006); Guieu et al. (2006). 

Classified as M5.5 bviWilking et al. (2001) 

" Classified as M5.75 bv lLuh rnan (2003) and as M6.5 bv lGuieu et~all (12001) 



TABLE 6 
M6 Sequence - Surface Gravities 





Age 


Prev. 


















Adopted 


Object 




Log(g) 


N3 


65 


64 


63 


62 


61 


59 


58 


Log(g) 


GY 5 


<1 Myr 


3.65 


3.80 


3.57 


3.30 


(4.14) 


(4.08) 


(4.39) 


(3.80) 


3.24 


3.75 


2MASS 2234+40AB 


~1 Myr 


4.00 


3.80 


3.54 


3.56 


(5.00) 


(4.81) 


4.19 


3.69 


3.60 


3.93 


CFHT Tau-7 


1.5 Myr 




3.50 


3.40 


3.44 


(4.00) 


(4.16) 


4.32 


3.80 


3.10 


3.72 


o-Ori 12 


3 Myr 


4.00 


3.70 


4.05 


3.59 


3.61 


3.74 


4.49 


3.78 


3.21 


3.78 


DENIS 1605-24 


5 Myr 


4.25 


3.70 


3.60 


3.49 


(3.60) 


4.47 


4.67 


3.94 


3.14 


3.82 


Gl 577BC 


70 Myr 


4.75 


4.43 


4.87 


(4.56) 


-(5.59)- 


(4.87) 


(5.09) 


(3.80) 


(4.07) 


4.68 


Gl 406 


>1 Gyr 


5.15-5.40 


5.67 


(5.61) 


(4.56) 


g go 


(4.85) 


5.60 




4^ 


5.19 



References. — Age and g ravity from ILuhman fc Riekd (ll999l):ILowrance et al.l ( 120051) ; IMcGovernI (120051 ) ; IMohantv et al.l 
( [2005D : IWilkine et al.l ( [20051) : [UaEallerol ( 120071) ; lAllers et al.l ( 120091) . 

Note. — See ^ 13.31 for explanation of annotations of the table entries. 



Rice et al. 



TABLE 7 
Velocity Results 







Previous 


Results 


Current 


Results 




Spectral 


RV 


V sini 


RV 


V sini 


Object 


Type 


(km s-l) 


(km s^^) 


(km s^-^) 


(km s-i) 


Upper Scorpius 


USco 66 AB 


M6 


-7.81,-4.4 


25.9, 27.5 


-5.8 


28 


DENIS 1605-24 


M6 






-5.2 


<7 


SCH 1612-20 


M6.5 






-6.8 


<9 


USco 100 


M7 


-7.2,-8.9 


43.7, 50 


-8.3 


50 


USco 130 


M7.5 


-3.91 


14, 15.2 


-7.1 


18 


SCH 1623-23 


M8 






-6.0 


53 


DENIS 1619-24 


M8 






-7.0 


47 


SCH 1622-19 


M8 






-10.2 


25 


TW Hydrae 


2MASS 1139-31 


M8 


11.6, 9.7 


25 


5.9 


30 


2MASS 1207-39AB 


M8 


11.2, 8.7 


13 


7.5 


18 


TWA 5B 


M8.5 


13.4 


16 


9.7 


27 


Other Young Objects 


2MASS 2234+40AB 


M6 


-10.6 




-13.1 


17 


CFHT Tau-7 


M6.5 






15.6 


20 


a Ori 12 


M6 


29.8, 37 




31.2 


18 


GY 5 


M6 


-6.39,-6.3 


16.5, 16.8 


-7.1 


20 


Gl 577BC 


M5.5 






-6.2 


43 


2MASS 0608-27 


M8.5 






22.6 


20 


I'ield Objects 


Gl 406 


M6 


19.1 


<2.9 


18.7 


<8 


LP 402-58 


M7 






-2.6 


17 


LP 412-31 


M8 


41.6, 44.7 


8,9 


41.8 


33 


2MASS 0140-1-27 


M9 


9.6, 8.2 


6.5 


8.6 


11 



References.^ — 'Tinney & Reid (ToM); Roid ct al.' ("2002); 'Reid' f2003]; 

[Moha ntv et all II2003I): [M ohantv & Basri (2003); MuzcroUc ct al. (2003); Kenyon et aL| 
Jmrat: Mohantv et aP \200S ): iKurosawa et al.l (|2006il : iReiners k Basril \200% : 
lAIIers et a l. (200^). 



